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The magnetic susceptibility of electrons confined to a 
spherical cavity or a circular billiard shows slow oscillations 
as a function of the number of electrons, which are a new 
manifestation of the Super Shell Structure found in the free 
energy of metal clusters. The relationship of the oscillations of 
the two different quantities is analyzed by means of semiclas- 
sical calculations, which are in quantitative agreement with 
quantal results. The oscillations should be observable for en- 
sembles of circular ballistic quantum dots and metal clusters. 

PACS numbers: 05.45.+b, 02.50.+s,03.65.-w,03.65.Ge 



I. INTRODUCTION 

The confinement of independent Fermions in 2 or 3 di- 
mensions (2D or 3D) leads to a bunching of the single 
particle levels, if the mean free path of the Fermions is 
large compared to the size of the system. This is known 
as shell structure (SS) and leads to oscillations of the 
total energy as function of particle number N around 
a smoothly changing background. The oscillating part, 
referred to as shell energy, has minima at the so called 
magic numbers, which have been known for a long time 
for nuclei. More recently they have been observed in 
the abundance spectra of alkali metal clusters (c.f. ||] 
and the original work cited therein), representing min- 
ima of the free energy 0. Later on, it has been found 
H , that the amplitude of these shell oscillations is mod- 
ulated by a slow oscillation. This so called Super Shell 
Structure (SSS) had been predicted theoretically ||,[|. 
In this paper we will show that the magnetic susceptibil- 
ity follows a similar SSS pattern. Using the Strutinsky's 
shell correction method and semiclassical Periodic 

Orbit Theory (POT), we will trace the SSS of the sus- 
ceptibility and free energy back to the same interference 
pattern between electrons on classical periodic orbits. 

The consequences of SS for the magnetic susceptibility 
of the confined electron gas have been discussed in refs. 
|5] |l3| and earlier references cited therein. The suscepti- 
bility of a confined ballistic 2D electron gas can be mea- 
sured for large ensemble of quantum dots on a AlGaAs - 
GaAs semi conductor hetero structure ||. For this type 
of experiments it is claimed [l3) that the shell os- 

cillations as a function of the electron number N are av- 
eraged out by the fluctuations of the size and the shape 
of the individual dots. The only Quantum Size Effect 



(QSE) expected to survive is a paramagnetic enhance- 
ment of the susceptibility, which changes smoothly with 
N. In this paper we will argue that such experiments 
should permit to resolve the slow oscillations reflecting 
the SSS. We will also discuss the experimental possibili- 
ties to detect SSS of the magnetic susceptibility of metal- 
lic clusters. 



II. SUSCEPTIBILITY OF ELECTRONS IN A 
SPHERICAL CAVITY 

Choosing z as the direction of the magnetic field H, 
the orbital part of the electronic Hamiltonian is |l5| ] 
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where Hq is the Hamiltonian at the zero magnetic field, 
consisting of kinetic energy and the confining potential, 
which is assumed to be circular or spherical. The opera- 
tor L z is the angular momentum projection on the z-axis 
and M the effective electron mass. We use the Larmor 
frequency uj = psH/h as the unit of the magnetic field 
H in order to stress the analogy to the case of a system 
rotating with the angular velocity — uj. 

Up to third order of perturbation theory in w, the ther- 
modynamical potential il(T, A, uj) as function of the tem- 
perature T and chemical potential A is 



n(r,A,w) = -TE I/ ln 
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where £„, m v and (. . .)„ are the energy, the angular mo- 
mentum projection and the expectation value with the 
unperturbed electron state v. In our units, the zero 
field susceptibility is a moment of inertia. For the grand 
canonical ensemble it reads 
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where n v = (1 + exp[(e 1/ — A)/T]) _1 are the Fermi occu- 
pation numbers and 9 ng is the moment of inertia of rigid 
rotation, p(r) being the particle density. 

The shell structure of 9 cr has been analyzed for nuclear 
rotation (where it is called cranking moment of inertia) . 
Following the concepts of Strutinsky's shell correction 
method fl, 9 cr and other quantities are divided into 
a smooth part and into an shell part (denoted by the 
subscript SH). The shell contribution to 9 cr , which we 
call 9sh, represents the total QSE, because the shell part 
of 9 r i g is negligible as compared to 9sh- The partition 
can either be done numerically starting from the quantal 
electron levels ||-|| or it can be based on semiclassical 
periodic orbit theory (POT). In the latter case the shell 
terms read [ p^| , p^|Jl^1 
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where £Ish denotes the value at zero field and a = 
in 2D or 3D, respectively. For the spherical cavity 
the circular billiard (l^Jl^l the orbits (3(t,p) are defined 
by the number t of the revolutions around the center and 
the number p of the corners, 
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where we have introduced the length of the orbit Lp, the 



wave number k = y2MX/h 2 and fp = 1 for diameters 
and 2 for planar orbits. The phases up, which are related 
to the Maslov indeces, are not important in our discussion 
and are given in [^,^3). The POT level densities of a 
spherical cavity and a circular billiard in a magnetic field 
have been studied in refs. |l^ , [l9[] . 

Since the energy to extract one electron from the con- 
fining potential is much higher than the temperature, it 
is important to use the canonical ensemble and define 
the susceptibility as the derivative 9 = —(d 2 F/duj 2 ) UJ= o 
of the free energy at fixed particle number N. The im- 
portance of the fixed electron number for the magnetic 



properties of 2D-structures has been pointed out previ- 
ously jl2],[l3| . We adopt the approximation valid for large 
N, calculating F(T, N, u) = Q(T, A, oj) + \N, where A is 
found from the condition 



dft(T, A, oj) /d\ = N 



(14) 



The zero field susceptibility is given by the expressions 
d) or (|) taken at \(N) fixed by eq. (juj) at w = 0. We 
solve this equation numerically both for the semiclassical 
and the quantum calculation. 



III. SPHERICAL CAVITY 

Fig. 1 shows 9sh(N,T) calculated by means of the 
numerical Strutinsky averaging procedure S from the 
quantal levels in a spherical cavity. The parameters are 
appropriate for sodium. As a unit we use the Landau 
diamagnetic susceptibility (LU) for the electron gas in the 
cavity, \9l\ — 0.2715MNrg. Since all contributions to 
the bulk susceptibility are of the same order of magnitude 
, the figure shows directly the enhancement due to the 
QSE. 




FIG. 1. Shell contribution 8s h of the magnetic sus- 
ceptibility of N electrons in a spherical cavity. The ef- 
fective mass is M — M e , the radius R = raA'' 1 / 3 with 
t s — 0.208 nm, the Fermi energy ef = 3.24 eV and the 
temperature T = 0.005ef = 170 K. The thin line shows the 
results without averaging over JV 1/3 and the fat ones after av- 
eraging with a Gaussian of width AjV 1,/3 = 0.75. Inset: solid 
and dashed lines show the quantum and semiclassical calcu- 
lations for the canonical ensemble, respectively. Dots present 
the grand canonical quantum result (fixed A = £f). 

The susceptibility 9sh oscillates with the period of the 
shells. In addition its amplitude is modulated with a slow 
oscillation, which is the "Super Shell Structure" (SSS), 
first noticed for the level density [[HJ. It was found ex- 
perimentally in the abundances of Na clusters |3| , which 
are determined by Fsh @ , showing also the SSS pattern 
(cf. Fig. 2). 

We have also evaluated 9$h by means of the POT 
sums (^|). As seen in the inset of Fig. 1, the quan- 
tal values of 9sh agree very well with ones obtained 
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from the semiclassical expression. Semiclassics permits 
a simple interpretation of the SSS. The shortest PO's en- 
closing magnetic flux are the triangle and square. The 
temperature factor, 1/ sinh(Trs/7i) w 2 exp(— Lp/Ax) , 
where At = h k/TM is the characteristic "tempera- 
ture length", damps the longer orbits. The beat pattern 
results from the superposition of the two leading terms, 
the triangle and the square. The basic oscillation has 
a period given by k(L& + L a )/2 and the beat oscillates 
with k(Ln — Fig. 3 shows a calculation that takes 

into account only the triangle and the square. Compar- 
ing with the full calculation, the influence of the longer 
orbits is seen. Although they make the peaks higher the 
beat pattern is not much changed up to N 1 ^ 3 ~ 16. The 
upper and lower envelops of the full calculation are not 
very different from the ones of the truncated calculation. 
This demonstrates that the beat pattern is basically gen- 
erated by the triangle and the square. For iV 1 / 3 > 16, the 
full calculation has a beat minimum where the truncated 
one has a maximum, indicating that the interference with 
the longer orbits becomes important. 
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FIG. 2. Shell contribution Fsh to the free energy of N 
electrons in a spherical cavity. The calculation and the con- 
ventions are identical with Fig. 1. 

IV. RELATION TO THE SHELL STRUCTURE IN 
THE FREE ENERGY 

The relation between the shell contributions to the zero 
field free energy, Fsh, and to the susceptibility, Osh, is 
understood by comparing the two POT sums (||). The 
terms are identical up to the factor 

= a(lf}Tff/h) 2 = 7C 2 M 2 T&ap 2 sin 2 (20) (15) 

in Osh, which suppresses the orbits lp — and gives 
the long orbits a higher weight. If only few orbits with 
similar values of Tp « T contribute, Osh ~ YFsh- The 
simple scaling is also expected to hold for shapes not too 
different from the sphere and can be used to relate the 
SS in the susceptibility and free energy. 

In the 3D case the diameter orbit is suppressed by a 
factor 1/ykR as compared to the planar orbits, because 



it has a lower degeneracy. This makes its contribution 
to F S h insignificant for large N. The two sums Fsh 
and 9sh become similar, because the leading terms are 
the triangle and square. In fact, for T > 0.02ef we 
find Osh ~ YFsh with p = 3 — 4 in ( |l5| ) (where £f 
is the Fermi Gas energy). For the lower temperature 
T = 0.005ef, shown in Figs. 1 and 2, we find a ratio T 
with p = 4 at the SSS maxima 

(N 3 = 12 and 17). The 
SSS minima are less pronounced for Osh than for Fsh- 
There the triangle and square cancel each other. The 
main contribution comes from longer orbits, which are 
much more important for Osh than for Fsh, preventing 
Osh from becoming as small as Fsh- This interpreta- 
tion is also supported by Fig. 3. The difference between 
the full and truncated calculations, which represents the 
contribution of the longer orbits, is just an up-shift of the 
upper envelop . 
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FIG. 3. Shell contribution 8sh of the magnetic suscepti- 
bility of N electrons in a spherical cavity. The thin line shows 
the same as fig. 1, the thick line corresponds to a calculation 
where only the triangular and square orbits in the POT sums 
([]) are taken into account. For N 1 ^ 3 < 6 the level density is 
not positive definite in this approximation and the solution of 
eq. ( |ll| ) is not possible. 

The SSS is clearly developed for the temperature of 
T = 0.005e F , shown in Fig. 1. For T = 0.0005e F 
we find very pronounced SS with an amplitude of a few 
thousand LU, but little SSS. The averaged susceptibility 
stays around 1500 LU with small peaks (400 LU high) 
at N 1 / 3, — 7 and 13. Due to a weak temperature damp- 
ing in the POT sum many orbits contribute to the sum, 
destroying the simple beat pattern. This is at variance 
with a distinct SSS pattern in Fsh persisting to T = 
The reason is the factor 1/T^ which suppresses the 
orbits with high p. For T = 0.05ef there is SSS, but the 
amplitude of the basic shell oscillations remains below 5 
LU for A^ < 100, becoming small compared with the LU 
for larger N. Hence, the intermediate temperature seems 
to be optimal for the observation of the SSS. 
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V. CANONICAL ENSEMBLE 

Experimentally, the average susceptibility of an ensem- 
ble of clusters with a distribution in N will be measured. 
Let us consider a Gaussian distribution in N 1 ^ 3 of width 
A(iW 3 ) — 0.75, which corresponds to about one os- 
cillation. As shown in Fig.l, averaging the susceptibil- 
ity damps out only the basic shell oscillations, whereas 
the SSS remains as a modulation of the strongly para- 
magnetic susceptibility. Hence, it is expected that the 
SSS can be observed with a moderate mass selection of 
the clusters that would not allow to resolve the basic SS 
(AN/N = 0.23 for N = 1000 in Fig. 1). 

The QSE of the susceptibility survives the averaging 
only for the canonical ensemble [g The inset of 

Fig. 1 compares the canonical with grand canonical en- 
semble (A = Ef)- Though the positions of the extrema 
are similar, the shape of peaks is rather different. For 
the grand canonical ensemble the negative and positive 
values are equally probable and, as the result, averag- 
ing with respect to N quenches 8sh- For the canonical 
ensemble, the positive values are more frequent and the 
QSE survives the averaging with respect to N. The pref- 
erence of 6 sb > for the canonical ensemble is evident 
from POT sum (^): The shell correction to the level den- 
sity gsH, which is given by setting the factor {...} = 1, 
oscillates in phase with 9sh, i- e. more N values cor- 
respond to a paramagnetic than to a diamagnetic sus- 
ceptibility and the average with respect to the particle 
number is positive. 

Let us discuss this correlation in more detail, because 
it is crucial for the measurability of the SSS. Fig. 3 illus- 
trates that the typical inverted parabolas appear already 
for the lowest orbits. It is sufficient to consider only one 
term in (^|), say the triangle. For the grand canonical 
ensemble, the wave number k = kp is constant (kp is the 
Fermi Gas value). The orbit length is La oc TV 1 / 3 . The 
susceptibility is proportional to sin(fci?LA) (for simplicity 
the phase is left away), which averages to zero. The 
particle number expectation value as given by (JTJ|) con- 
tains a term proportional to — cos(kpL&) that makes it 
oscillating around N. Thus, for the canonical ensemble, 
the wave number k = kp+Sk cannot be constant. It must 
contain an oscillating term 5k in order to satisfy eq. (14). 



In order to understand qualitatively its consequences for 
the susceptibility we use an argument from refs. |l2| , |l3| . 
Retaining only the linear order of Sk, eq. (M) gives 
8k oc cos{kpL/\). The susceptibility is proportional to 
sin((kp + 5k)L&) ~ sm(kpL/\) + 5kL& cos(kpL/\). 
The second term is proportional to cos(kpL&) 2 which 
averages to 1/2. Hence, the total averaged susceptibility 
is positive. Since the interference between the triangular 
and quadratic orbits in both 9sh and gsn is about the 
same, the SSS modulates paramagnetic term also after 
averaging. The oscillations of Sk narrow the minima and 
broaden the maxima of the susceptibility. In our calcu- 
lations, Sk is treated exactly. As seen in Figs. 1 and 3, 



due to the higher orders in Sk the minima are narrowed 
to cusps and the maxima take the shape of parabolas. 



VI. REAL CLUSTERS 

Real clusters deviate from the perfect spherical cav- 
ity: The surface has a finite thickness of the order of the 
screening length. The discrete ionic back ground implies 
a certain surface roughness of the order of the interatomic 
distance. There may be impurities or other imperfections 
distributed over the volume. 

The SS in a spherical potential with a realistic surface 
thickness (the one of sodium) has been studied in ref. 
@. The SSS in the binding energies is clearly developed. 
The beat minima are somwhat shifted as compared to 
the cavity. The shift has been be traced back to small 
changes of the action of the POT due to the modified "re- 
flection" by the finite potential at the surface. Since the 
expressions (JgJ) hold also for the more realistic potentials 
if the appropriate action S is inserted, a similar shift of 
the SSS pattern can be infered for the susceptibility. 

The consequences of the surface roughness for the shell 
structure of the ground state energy have been studied 
in ref. pi| . The energy is given by the POT sum (||), 
where each term contains an additional damping factor 
X p '. Assuming that the rough surface is randomly dis- 
placed relative to the ideal one, the damping factor is 



X 1 



exp {—2p(ak sin < 



(16) 



where a Gaussian displacement distribution with the 
width a ~ rs is used. It has a simple interpretation. 
The number of reflection on the surface is p. For each re- 
flection the rough surface scatters a certain fraction of the 
particles away from the POT. The arguments of ref. pl| ] 
can be immediately taken over to the zero field suscep- 
tibility. The damping factor arises from the reflections 
on the irregular surface. A weak magnetic field does not 
change the reflections on the surface [] and the damping 
factor (16) appears also in the susceptibility. Hence, the 
long orbits (large number of reflections p) are strongly 
suppressed by the surface roughness. In ref. [|IJ it is 
shown that a roughness a — 0.2rs — 0.38//cf reduces 
the amplitude of the SS to about 1/2 of the one of the 
ideal cavity. The SSS is found to be nearly the same 
as for the ideal cavity. For the susceptibility one expects 
that the surface roughness strongly reduces the contribu- 
tions of the long orbits, such that the SSS pattern of Fig. 
1 approaches the one of Fig. 3 with a reduced amplitude. 



The difference beween the case with a weak and without 
a field is the small curvature of the trajectory between the 
reflection points which barely changes the angles of the tra- 
jectory with the surface. 
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In particular, for very low temperature the surface rough- 
ness is expected to enhance the SSS, because it efficiently 
damps the long orbits. 

Impurities or other imperfections that are homoge- 
neously distributed over the volume will also scatter the 
particles away from the POT. The situation is analogous 
to the propagation of a wave in an absorbing medium, 
which has been considered in ref. The scattering 

results in a damping factor of the form exp (— Lp/Aj), 
where 1/Aj measures the amount of scattering per unit 
length and corresponds to the mean free path of the elec- 
trons due to the imperfections. This kind of damping is 
equivalent with an increase of the temperature, because 
the temperature damping factor in 1/2 sinh (Tt^/Ti) 
exp (— hp I At). The consequences of a temperature in- 
crease are discussed above. 



VII. MEASUREMENTS OF THE 
SUSCEPTIBILITY OF METAL CLUSTERS 

In Ref. |H the susceptibility has been measured for 
three ensembles of clusters with the average size of 
TV 1 / 3 »: 4, 2 and 1.5 and a large spread in size. A param- 
agnetic enhancement of 5, 2 and 1.5 is found, respectively. 
The decrease of the enhancement with N, which ref. H] 
mentions as an unexplained phenomenon, can be seen 
in Fig. 1. To identify the SSS, experiments with more 
points in N and a mass resolution better than 30% are 
needed. An alternative experiment would be the mea- 
surement of the deflection of a cold cluster beam in an 
inhomogeneous magnetic field, which provides directly 
the susceptibility. To reach the necessary sensitivity of 
such a Stern-Gerlach apparatus seems to be possible p2[ . 
It is favourable for this kind of experiment that only a 
moderate mass selection is needed to observe SSS, what 
permits larger intensities. 

Measuring the susceptibility could shed new light on 
the electronic structure of small metal clusters. Solid 
icosahedral shapes have been suggested for sodium clus- 
ters with T ~ 200AT Our sphere model should be 
a rough first approximation for N < 1000. In this range 
the shell energies Esh for the spherical and icosahedral 
cavity are similar (24) and, as discussed above, the same 
can be expected for the susceptibilities. A more pro- 
nounced paramagnetic SSS is expected if shape of the 
cold clusters comes close to a rough sphere. The picture 
changes, if the clusters were liquid or would keep keep 
the same shape as in the liquid state when freezing at 
low T. Then magic clusters are spherical and strongly 
diamagnetic, whereas the non magical ones are deformed 
and weakly diamagnetic j2(J. Thus, the averaged suscep- 
tibility would be diamagnetic, still showing a SSS. 



VIII. HALF SPHERE 

Fig. 4 shows the susceptibility for a half sphere with 
R = r s (N/2)^ 3 . The results are very similar to the full 
sphere. The main difference is a shift of the basic shell os- 
cillations by half a period. The averaged susceptibility is 
almost indistinguishable from Fig. 1. Cluster deposited 
on an insulating surface may take shapes close to half 
spheres M-BTj. 




FIG. 4. Shell contribution of the magnetic susceptibility N 
electrons in a half spherical cavity. The line conventions and 
calculation (except R — rsiN/2) 1 ^ 3 ) are identical with Fig. 
1. 



IX. CIRCULAR QUANTUM DOTS 

The susceptibility of the 2D electron gas confined to a 
circular potential well with R = 500 nm is shown in Fig. 
5. An effective mass of M = 0.067M e pq ], appropriate 
for GaAs, is assumed. The LU is \6 L \ = MR 2 /?,. The 
2D case is similar to the 3D - case, the main difference 
consisting in an increase of the QSE with N (The val- 
ues in Fig. 5 are divided by N). At the temperature 
T = 5h 2 /2MR 2 « 0.13 AT a distinct SSS is seen. For 
T = 0.5h 2 /2MR 2 the averaged susceptibility ~ 4/N LU 
showing shallow oscillations with an amplitude of ~ 1 /N 
LU, which loosely correlate with the SSS in Fig. 5. A SSS 
beat pattern for a circular dot has been first calculated 
|l3[| . However, it is argued there and also in [fj"l 12 



that when averaging over the ensemble of dots used as 
experimental probe, uncertainties in the shape and size 
will completely wipe out the shell structure, the only re- 
maining QSE being a paramagnetic enhancement that 
varies smoothly with N . In contrast, Fig. 5 shows that 
averaging with a Gaussian of width IAN 1 / 2 = 1.6 ( cor- 
responding to a 10% spread in N for N = 1000) destroys 
only the basic SS, whereas the SSS remains visible. 

Nowadays it is possible to manufacture probes with a 
large number of circular quantum dots, specifying the 
radius and the gate voltage with a 5% accuracy [p8[ . 
Changing the number of electrons in the dots by means of 
the gate voltage seems to be a possibility to measure the 



5 



SSS of the susceptibility. The imperfections in manufac- 
turing circular dots will have similar effects as discussed 
above for the metal clusters. The discussion of the 3D 
case can directly be applied to the 2D case. Thus, the 
SSS pattern is expected to survive if the surface rough- 
ness a < 0.4/kp. 




FIG. 5. Shell contribution 8sh to the magnetic suscepti- 
bility of N electrons in a circular well. The effective mass is 
M — 0.067M e , the well radius R = 500 nm and the tempera- 
ture T = 5h 2 /2MR 2 » 0.13K . The line conventions are the 
same as in Fig. 1. Averaging over N 1 ^ 2 is carried out with 
a Gaussian of width AN 1/2 = 1.6. The non averaged values 
are divided by a factor of 10. 




FIG. 



Shell contribution Fsh to the free energy of N 



electrons in a circular well in units of ef = Nh /MR . The 
calculation and the conventions are identical with Fig. 3 

Discussions with M. Brack and F.A.Ivanyuk and finan- 
cial support by INTAS (93-0151) are gratefully acknowl- 
edged. 



Fig. 6 demonstrates that, at variance with the 3D 
case, Fsh{N) significantly differs from 0sh(N). For 2D 
case, the diameter orbit is not suppressed in Fsh- The 
interference of the diameter, triangle and square results 
in slower basic and faster beat oscillations as compared 
to 6sh, for which the diameter is missing. Averaging 
with respect to iV 1 / 2 filters out the slow oscillation due 
to the interference between triangle and square, which 
mainly modulates Fsh(N). The SS of the free energy 
should show up as a modulation of the capacitance of 
the dot. It is given by d 2 F/dN 2 , which also determines 
the abundances of heavy clusters . 



X. CONCLUSIONS 

The susceptibility of electrons confined in two or three 
dimensions by a spherical potential oscillates as func- 
tion of their number. This shell structure is modulated 
by slow oscillations, the Super Shell Structure, which 
only develops at sufficiently high temperatures. Mea- 
surements that average out the shell structure may still 
reveal the Super Shell Structure. The free energy of elec- 
trons confined in three dimensions, shows the analogous 
Super Shell pattern, which is observed in the abundances 
of metal clusters. However, for the two dimensional po- 
tential the shell structure of the free energy differs con- 
siderably from the one of the susceptibility. 
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